HNCO in massive galactic dense cores

A&A manuscript no.

(will be inserted by hand later)

Your thesaurus codes are:
08(08.06.2; 09.03.1; 09.13.2; 13.19.3)

ASTRONOMY

AND
ASTROPHYSICS

HNCO in massive galactic dense cores *

I. Zinchenko^ 2, c. Henkel^, and R.Q. Mao^^

^ Institute of Applied Physics of the Russian Academy of Sciences, 46 Ulyanov str., 603600 Nizhny Novgorod, Russia
^ Helsinki University Observatory, Tahtitominmaki, RO.Box 14, FIN-00014 University of Helsinki, Finland

Max-Planck-Institut fiir Radioastronomie, Auf dem Hiigel 69, Bonn, Germany
* Purple Mountain Observatory, 210008 Nanjing, China

Received date ; accepted date

Abstract. We surveyed 81 dense molecular cores associated
with regions of massive star formation and Sgr A in the
•^K-iKi = 5o5 — 4o4 and lOgio — 9o9 lines of HNCO. Line
emission was detected towards 57 objects. Selected subsam-
ples were also observed in the Iqi Oqo, 4o4 — 3o3, 7o7 — 6o6,
15oi5 - 14oi4, I6016 - 15oi5 and 2I021 - 2O020 lines, cover-
ing a frequency range from 22 to 461 GHz. HNCO lines from
the K-i = 2,3 ladders were detected in several sources. To-
wards Orion-KL, K-i = 5 transitions with upper state ener-
gies E^/k ~ 1100 and 1300 K could be observed.

Five HNCO cores were mapped. The sources remain spa-
tially unresolved at 220 and 461 GHz (lOoio — 9o9 and 21oio —
2O020 transitions) with beam sizes of 24" and 18", respectively.

The detection of hyperfine structure in the Iqi — Oqo transi-
tion is consistent with optically thin emission under conditions
of Local Thermodynamic Equilibrium (LTE). This is corrobo-
rated by a rotational diagram analysis of Orion-KL that indi-
cates optically thin line emission also for transitions between
higher excited states. At the same time a tentative detection of
interstellar HN^CO (the 10o,io - 80,9 hne at 220 GHz towai'd
G 310.12-0.20) suggests optically thick emission from some
rotational transitions.

Typical HNCO abundances relative to H2 as derived from
a population diagram analysis are ~ 10~^. The rotational tem-
peratures reach ^ 500 K. The gas densities in regions of HNCO
K^i = emission should be n ^ 10^ cm^-^ and in regions of
K_i > emission about an order of magnitude higher even
for radiative excitation.

HNCO abundances are found to be enhanced in high-
velocity gas. HNCO integrated line intensities correlate well
with those of thermal SiO emission. This indicates a spatial

Send offprint requests to: I. Zinchenko

* Based on the observations collected at the European Southern Ob-
servatory, La Silla, Chile and on observations with the Heinrich-Hertz-
Telescope (HHT). The HHT is operated by the Submillimeter Tele-
scope Observatory on behalf of Steward Observatory and the MPI fiir
Radioastronomie.

** Tables 1, 2, 5, 6 are also available in electronic form and Ta-
bles 7-14 are only available in electronic form at the CDS via anony-
mous ftp to cdsarc.u-st rasbg.fr (130.79.128.5) or via h ttp://cdsweb.u-
strasbg.fr/Abstract.html
Correspondence to: zin@appl.sci-nnov.ru

coexistence of the two species and may hint at a common pro-
duction mechanism, presumably based on shock chemistry.

Key words: Stars; formation - ISM: clouds - ISM: molecules
- Radio lines: interstellar

1. Introduction

Systematic studies of dense molecular cores in regions of high
mass star formation (HMSF) are of great importance for our
general understanding of star formation. In comparison with
low mass star formation regions, so far only a few rather arbi-
trarily selected cores associated with HMSF have been investi-
gated in some detail.

In recent years we performed extensive surveys of dense
cores in regions of high mass star formation, mainly in CS
(Zinchenko et al. 1995, 1998). We used water masers as sign-
posts of high mass star formation. Both outer and inner Galaxy
were covered by these surveys (I « 120° — 210° and 260° —
308°). The innermost part of the Galaxy (I « 308° - 360°)
was observed in a similar way by Juvela (1996). In addition,
sources as sociated with wa ter masers were surveyed in ther-
mal SiO ( ^-[arju et al. 1998 ) which is supposed to be a good
indicator of shocks in molecular clouds. From these obser-
vations we derived basic physical parameters of the cores

and constructed their statistical distributions (Zinchenko 1995
^inchenko et al. 1998 ). In order to investigate a range of core

densities, observations of lines with different excitation condi-
tions are needed. One of the interesting candidates is the HNCO
(isocyanic acid) molecule.

HNCO was first detected by Snyder & Buhl (1972) in
Sgr B2. Subsequent studies have concentrated mostly on the
Galactic center region where the HNCO emission was found
to be particularly strong (e.g., Churchwell et al. 1986, Wilson
et al. 1996, Lindqvist et al. 1995, Kuan & Snyder 1996, Dah-
men et al. 1997, Sato et al. 1997). A survey of HNCO emis-
sion throughout the Galaxy was made by Jackson et al. (1984)
in the Jk_iKi — 5o5 — 4o4 and 4o4 — 3o3 transitions with
the 11m NRAO telescope. Seven (from 18) clouds including
Orion KL were detected at rather low levels of intensity (typi-
cally - 0.2 K on a TX scale). Churchwell et al. (1986) obtained

2

I. Zinchenko et al.: HNCO in massive galactic dense cores

strict upper limits on HNCO Iqi — Oqo and 2o2 — loi emission
towards about 20 galactic sources with the 36.6 m Haystack
antenna.

HNCO is a slightly asymmetric rotor Its levels may be des-
ignated as Jk_iKi where J is the total angular momentum and
K-i, Ki are quantum numbers corresponding to the projection
of J on the symmetry axis for the limiting cases of prolate and
oblate symmetric top, respectively (e.g. Townes & Schawlow
1975). The structure of the HNCO energy levels can be rep-
resented as a set of "ladders" with different K^i values, like
for a symmetric top. However, due to the asymmetry of the
molecule radiative transitions between different K^i ladders
(5-type transitions) are allowed and, moreover, they are very
fast. The corresponding component of the dipole moment is
similar to its component for transitions inside the K^i ladders
(a-type transitions). Churchwell et al. (1986) found that as a re-
sult the HNCO excitation is governed mostly by radiative rather
than collisional processes (at least in Sgr B2).

On the basis of their estimates of source parameters Jackson
et al. (1984) concluded that HNCO is a potentially valuable
probe of the densest regions (n ^ 10^ cm^"^) of molecular
clouds. It was shown also that HNCO is rather sensitive to far
infrared (FIR) radiation fields due to the fact that the lowest
levels of the K^i = 0, 1 and A'_i = 1,2 ladders are separated
by energies corresponding to FIR wavelengths (330 /im and
110 fim, respectively).

From this consideration it is clear that multitransitional data
are needed to understand HNCO excitation and to derive the
source properties. Bearing this in mind we undertook a survey
of HNCO emission in various rotational lines, also trying to
detect emission from higher excited ladders (K^i > 0).
Five cores were mapped in HNCO to estimate the extent of the
emission.

Several other species were observed simultaneously with
HNCO. The most prominent are C^*0 and SO. In the following
we thus also compare HNCO with C^'^O.

2. Observations

2.1. Source list

For this study we observed those dense cores showing par-
ticularly strong CS emission (Tmb > 3 K) in the surveys of
Zinchenko et al. (1995, 1998) and Juvela (1996). Several strong
SiO {v = 0) sources detected by Harju et al. (1998) are also
included in our sample. Sources observed at the SEST and at
Onsala are presented in Tables |l], |2| Sources also observed at
Effelsberg or at the HHT are marked in both tables.

We designate most sources according to their galactic co-
ordinates. Exceptions are Orion KL and Sgr A. For Sgr A we
use the position observed by Jackson et al. (1984) for com-
parison (known as the M-0. 13-0.08 cloud, see Lindqvist et al.
1995). Common identifications with some well known objects
are given in the last column.

Table 1. Source list for SEST observations.

Name

a(1950)

5(1950) Remai-ks
(°) (')(")

o zoi.o4— z.uy

AO QA OQ O
Uo JU Zj.Z

A Q AQ Q 1
— 4J UJ Jl

OA/1 OO 1 1 AQ

no c/1 m A
Uo j4 jy.U

/I O CQ QA
— 4Z jj JU

Dr^"\l7 Q/1

KL.W J4

tj Zoj. 14+1.43

no C7 1A 1
Uo J / JO. J

Al QQ Q

— 4J JJ Jo

KL,W JO

^ o/;'7 c\A 1 A/;
(j Zo /.y4— l.Uo

no C7 O 1 T

Uo J / Z 1 . /

AH 1 A f\A

—4/ ly U4

KL.W Jo

U Z0O.4Z— U.OJ

nn nn i o i
UV UU IZ.l

AH QO m

—4/ JZ U/

1 A. 11/1
O ZoV. Id— 1. 14

nn n 1 c i a
UV Ul J 1.0

AQ 1 A /I Q

— 4o 10 4J

O Z/U.Zo+U.oi

nn ^ A con
Uy 14 D5.U

An A A nn
—4 / 44 UU

KLW 41

Lj Zoj.Zd— U.Uj

1 A on 1A o
lU Zy jD.o

CT /1A /1A

— J / 40 4U

U ZoD.ZUH-U. 1 /

1 A Q A Q /I O

lU Jo J4.0

CO AQ OO
— JO UJ ZZ

^ on 1 T7 ATI

U Z91.Z/— U. / 1

I 1 nn /ion

I I UV 4Z.U

A 1 A 1 CC

— ol Ul JJ

OA 1 r\ A'2

11 1 o c/1 n
11 IZ j4.U

AA CO C?

— dU jZ J /

"NT/^/^ QAAQ

INOL. JOUJ

O Z94.y/ — 1./3

11 1 A C 1 A

11 Jo jl.O

AQ 1 o nn
— dJ IZ UV

QAA m 1 1 1 /I

O JUU.y / + 1.14

1 o 10 An o
IZ JZ UU.Z

A 1 OQ A A

— ol ZJ 44

Tll^WT AC

KUW oD

O jUI. IZ— U.ZU

1 O 10 Q 1 "Q

IZ jZ ji.j

AO A A Q

— OZ 44 Jo

ir\z on 1 A oi
U JUj.ZU+U.ZI

1 Q m CA A

IJ U/

AO 1 O CA

— oZ lo DU

KCW /4

^ ir\^ 1 A oi

1 1 nn o 1 o
IJ UV Zl.Z

AO 1 O AO

— oZ lo UZ

^Ao on A o^

1 1 1 H OT O

IJ IJ Z/.Z

AO /10 CA

— dZ 4Z Do

QAO AA 1 O AO

O JUo.UU+Z.UZ

1 Q on O /I T

1 J Zy Z4.J

AA 1 1 OO

— dU 11 ZZ

QAO AO 1 A 1 O

u juo.yz+u. iz

1 Q 1C\ 1 A A
IJ jy J4.4

A 1 CQ /I C

— ol DJ 4D

KCW ly

ifWi no \ C\ AQ

1 Q /n 1 o c
IJ 4/ IZ.j

A 1 1 A CO

— ol Ly Do

iy c nn a ao
(j Mb. 1 /— U.UZ

1/1 /II 1 n /I
14 41 1U.4

C A Q C Q A

— jy JJ JU

1/1/11 QA A

14 41 Jo. 4

CA QA CQ

— jy JO DJ

1^Q A^ 1 A AA

1/1 /in c 1 A
14 4y J 1.5*

CO CA /1A

— Do DO 4U

QOQ HA A OC

U MS. /4— U.Z!)

1 C OT /I A O

1 J Z/ 4b'.o

CA on 1 c
— Do ZU ID

I'^A on 1 A 1 o
Lj jZ4.ZU+U. IZ

1 c on Al o
ij zy Ul.Z

CC /I A 1 O

— DD 40 IZ

'20A AH 1 A nr\

1 C QA OO O

Ij jy Zo.Z

CQ CO ni
— DJ Do Ul

QOA A/1 1 A Al

O JZd.d4+U.d1

1 C /I A /I O A

1 J 4U 4Z.D

CQ CA OA

— DJ Do Zy

Qoo in \ r\ Al

1 C CA ICQ

Ij jU id. J

CQ AO /I A

— D J UZ 4d

HOC 1 1 A Al

1 c c 1 CA n

Ij jl jy.V

CO Q/1 0/1

— DZ J4 Z4

/^QOOO/l AC/1

U JZ0.Z4— U.j4

1 C C/1 A/1 A

1 J j4 U4.9

CQ cn nn
— D J DU UV

"jon Al A on

1 c CA /in 1
1 J Do 4U. 1

CQ A/1 no
— DJ U4 Uo

Z"' QQAAC A 1A

1 A nA AQ A

lo Uo UJ.4

CI An Q A

— Dl 4/ JU

QQA OO A QA

U JoU.oo— U.Jo

1 A AA Qn n
lo Uo JU.U

C 1 CO 1/1

— Dl Do 14

n 1 OO A 1 Q

1 A m QA n
Id U/ Jo.U

CI QQ /I A

— D 1 J J 4U

1'20 0'2 A CC

1 A 1 A OQ n

Id Id ZJ. /

cn /I c /I c
— DU 4D 4D

D fWJ 1 AA

KCW lUo

111 11 r\ Al

1 A IT 1 O A

Id 1 / IZ.o

cn OO 1 o
— DU Zo lo

U jjJ.oU— u.zz

1 A 1 O 0/1 C

Id lo Z4.D

A A c A no
— 4V DV Uo

'!i'^'7 /in A /in

1 A Q C AO 1

lo JJ Uo.l

AH OO OQ

—4/ ZZ ZJ

1Ar\ AA A oc

1 A A A Q A A

lo 44 Jo. 4

/I C 1 A OA

— 4D lo Zo

1/1 c A 1 1 1 on
O j4j.U1 + 1.oU

1 A CQ 1 O O

Id jJ io.o

/in AA QA

— 4U UV Jo

a/iQ 1 n riA
O j4j.1Z— U.Uo

1 A C/1 /1 Q

10 J4 4Z.O

Al An /1Q
— 4Z 4/ 4V

T/l C CI 1 A

O j4j.j1+U.jj

1 / UU jJ.O

/in /in no

— 4U 4U UZ

n nn n 9^

O j^j.UU — U.Zj

1 7 ni An 7
1 / Ul ^u. /

Al n7
—"41 ZD U /

G 345.41-0.94

17 06 02.3

-41 31 44

G 348.55-0.97

17 15 53.1

-39 00 57

G 350.10+0.09

17 16 01.0

-37 07 30

G 348.73-1.04

17 16 39.7

-38 54 17

RCW 122

G 351.41+0.64

17 17 32.5

-35 44 13

G 351.58-0.36

17 22 04.2

-36 10 11

G 351.78-0.54

17 23 20.9

-36 06 53

G 353.41-0.36

17 27 06.5

-34 39 41

G 359.97-0.46

17 44 10.4

-29 11 03

Orion KL" ''

05 32 47.0

-05 24 23

G 173.48+2.45'=

05 35 51.3

35 44 16

S231

G 192.60-0.05''

06 09 58.2

18 00 17

S255

Sgr A"

17 42 28.0

-29 04 01

" also observed in Effelsberg, *
served in Onsala.

also observed with HHT, also

I. Zinchenko et al.: HNCO in massive galactic dense cores
Table 3. Observing parameters (Ai/ is the spectral resolution).

3

Molecule Transition" Frequency Telescope Date HPBW r/mb Tjys

(MHz) (") (K) (kHz)

HNCO

l()i — Ooo

21981.460

Eff. 100m

1998

40

0.3

50-100

12.5

4o4 — 3o3

87925.252

OSO 20m

1997

40^=

0.60"

210-290

250

505 — 4o4

109905.758

OSO 20m

1997

35"

0.52"

300-450

250

5o5 — 4o4

109905.758

SEST 15m

1997

47c

0.71"

200-270

86

7o7 — 6o6

153865.080

SEST 15m

1997

33"

0.64"

150-180

86

lOo.io — 9o,9

219798.320

SEST 15m

1997

24"

0.52"

190-360

86

15o,15 - 14o,14

329664.535

HHT 10m

1999

25

0.50

900-6000

480

16o,15 — 15o,14

351633.457

HHT 10m

1999

24

0.50

700-1000

480

21o,21 — 20o,20

461450.670

HHT 10m

1999

18

0.38

900-2000

480

1-0

109782.160

OSO 20m

1997

35"

0.52"

300-450

1000

2-1

219560.319

SEST 15m

1997

24"

0.52"

190-360

1400

"The frequencies and spectral resolutions for the observed A"_i > transitions are presented in Table ^
'The system temperatures are given on a scale.

"Beam sizes and main beam efficiencies are obtained by interpolating the data from the SEST manual (for SEST) and those provided by
L.E.B. Johansson (for OSO) at nearby frequencies.

Table 2. Source list for Onsala observations.

Name q(1950) (5(1950) Remarks

(") C") n n (') (")

G

121.30+0.66

00

33

53,

.3

63

12

32

RNO IB

G

123.07-6.31

00

49

29,

.2

56

17

36

NGC 281

G

133.69+1.22

02

21

40,

.8

61

53

26

W3 (1)

G

133.95+1.07*

02

23

17,

.3

61

38

58

W3 (OH)

G

170.66-0.27

05

16

53,

.6

36

34

21

IRAS05 168+3634

G

173.17+2.35

05

34

35,

.9

35

56

57

IRAS05345+3556

G

173.48+2.45'''"

05

35

51

.3

35

44

16

S231

G

173.72+2.70

05

37

31,

.8

35

40

18

S235

G

188.95+0.89

06

05

53,

.7

21

39

09

S 247

G

34.26+0.15

18

50

46,

.4

01

11

10

IRAS 18507+0 110

G

40.50+2.54

18

53

47,

.0

07

49

26

S76E

G

43.17+O.Or'''

19

07

49,

.8

09

01

17

W49N

G

49.49-0.39"'*'

19

21

26,

.2

14

24

44

W51 M

G

60.89-0.13

19

44

14,

.0

24

28

10

S87

G

61.48+0.10

19

44

42,

.0

25

05

30

S88B

G

70.29+1.60

19

59

50,

.0

33

24

17

K3-50

G

69.54-0.98''

20

08

09

.9

31

22

42

ONI

G

77.47+1.77

20

18

50,

.0

39

28

45

JC20188+3928

G

75.78+0.34

20

19

51

.8

37

17

01

0N2N

G

81.87+0.78''

20

36

50,

.5

42

27

01

W75 N

G

81.72+0.57"'''

20

37

13,

.7

42

12

11

W75 (OH)

G

81.77+0.60

20

37

16

.6

42

15

15

W75 S3

G

92.67+3.07

21

07

46,

.7

52

10

23

J21078+5211

G

99.98+4.17

21

39

10,

.3

58

02

29

IRAS21391+5802

G

108.76-0.95

22

56

38,

.4

58

31

04

JC22566+5830

G

108.76-0.98

22

56

45,

.2

58

29

10

S152(OH)

G

111.53+0.76"'''

23

11

36,

.1

61

10

30

S 158

" also observed in Effelsberg, '' also observed with HHT, " also ob-
served with SEST.

2.2. Observational procedures

The most important parameters of our SEST- 15m, OSO-20m,
Effelsberg 100-m and HHT measurements are summarized in
Tables ||, Further details are given belo'w for each instrument.

Table 4. Frequencies and spectral resolutions for the observed
K^i > transitions. For other observing parameters see Ta-
ble |

Transition Frequency Ai/

(MHz)

(kHz)

102,9 -

-92,8

219733,

.850

1400

102,8 -

-92,7

219737,

.193

1400

103,8 -

-93,7

219656

.710

1400

103,7 -

-93,6

219656

.710

1400

104,7 -

-94,6

219547,

.082

1400

104,6 -

-94,5

219547,

.095

1400

105,6 -

-95,5

219392,

.412

1400

105,5 -

-95,4

219392,

.412

1400

152,14

- 142,13

329573,

.46

480

152,13

- 142,12

329585,

.09

480

212,20

- 202,19

461336,

.93

480

212,19

- 202,18

461368,

.88

480

213,18

- 203,17

461182,

.51

480

213,19

- 203,18

461182,

.45

480

214,17

— 204,16

460950

.89

480

214,18

- 204,17

460950

.89

480

215,16

- 205,15

460625,

.75

480

215,17

— 205,16

460625,

.75

480

2.2.1. SEST observations

The observations ■were performed ■with SIS receivers in a
single-sideband (SSB) mode using dual beam s^witching ■with
a beam thro^w of ~ 12'. At 220 GHz ■we used 2 acousto-optical
spectrometers in parallel: (1) a 2000 channel high-resolution
spectrometer (HRS) ■with 86 MHz bandwidth, 43 kHz channel
separation and 80 kHz resolution and (2) a 1440 channel low-
resolution spectrometer (LRl) with a 1000 MHz total band-
width, 0.7 MHz channel separation and 1.4 MHz spectral reso-
lution. The LRl band was centered on the HNCO lOo.io — Qo.g

4

I. Zinchenko et al.: HNCO in massive galactic dense cores

219.4 219.5 219.6 219.7 219.8 219.9 220.0 220.1 220.2
Rest frequency (GHz)

Fig. 1. A SEST low resolution spectrum. For identified fea-
tures, molecular species and transitions are given

transition. However, it covered some other HNCO transitions
too (see Tableg) as well as C^^O (2-1), SO (65 - 64) and other
lines (Fig. |l] shows a typical spectrum).

The 110 and 154 GHz observations were performed simul-
taneously; the spectra were recorded by the HRS which band
was split into two equal parts. The 220 GHz HRS spectra were
smoothed to 170 kHz resolution and the 110 and 154 mm spec-
tra were smoothed to 86 kHz resolution. Pointing was checked
periodically by observations of nearby SiO masers; the point-
ing accuracy was ^ 5".

The standard chopper- wheel technique was used for the cal-
ibration. We express the results in units of main beam bright-
ness temperature (Tmb) assuming the main beam efficiencies
(Vmh) as given in Table |[ The temperature scale was checked
by observations of Orion KL.

In most sources only one position was observed, corre-
sponding typically to the peak of the CS emission. In addition,
G 270.26h-0.83 and G 301.12-0.20 were mapped with a spac-
ing of 10".

2.2.2. Onsala observations

At Onsala, the 110 GHz observing procedure was very similar
to that at the SEST. The observations were also performed in
a dual beam switching mode with a beam throw of 1 1.'5. The
front-end was a SIS receiver tuned to SSB operation. As back-
end we used 2 filter spectrometers in parallel: a 256 channel
filterbank with 250 kHz resolution and a 5 12 channel filterbank
with 1 MHz resolution. The calibration procedure was the same
as at the SEST. The pointing accuracy checked by observations
of nearby SiO masers was ^ 5". The strongest HNCO source
from the Onsala sample, W51M, was mapped with 40" spac-
ing.

2.2.3. Effelsberg observations

The 22 GHz observations in Effelsberg were performed with
a K-band maser amplifier using position switching. The offset
positions were displaced by 10'-15' symmetrically in azimuth.
Pointing was checked periodically by observations of nearby
continuum sources; the pointing accuracy was ^ 10". The in-
tegration time per position was a few hours.

The main beam temperature scale was checked by observa-
tions of nearby continuum calibration sources, NGC 7027 and
W3(OH); for Sgr A we used Sgr B2. The fluxes for the first two
sources were taken from Ott et al. (1994). The Sgr B2 flux at
1.3 cm was taken from Martm-Pintado et al. (1990).

2.2.4. HHT observations

To observe the HNCO J = 21-20 lines at 461 GHz we have
used the Heinrich Hertz Telescope (HHT) on Mt. Graham
(Baars & Martin 1996) during Feb. 1999 with a beamwidth
of 18". Spectra were taken employing an SIS receiver with
backends consisting of two acousto optical spectrometers with
2048 channels each, channel spacing —480 and —120 kHz, fre-
quency resolution —930 and 230 kHz, and total bandwidths
of —1 GHz and 250 MHz, respectively. Receiver temperatures
were — 150K, system temperatures were — lOOOK on a
scale. The receiver was sensitive to both sidebands. Any imbal-
ance in the gains of the lower and upper sideband would thus
lead to calibration errors. To account for this, we have observed
the CO J = 4-3 line of Orion-KL with the same receiver tuning
setup and obtain — 70 K, in good agreement with Schulz et
al. (1995).

HNCO J = 16 - 15 (351.63346 GHz) and J = 15 - 14
(329.66454 GHz) line emission was observed with a dual
channel SIS receiver in early April 1999 at the HHT. The
beamwidth was 22", receiver temperatures were 135 K; sys-
tem temperatures were —700 K on a scale. The receivers
were also sensitive to both sidebands. We have used pub-
lished spectra from Orion-KL and IRCh-10216 as calibrators
(proesbeck et al. 1994 ISchilke et al. 19971).

All results displayed are given in units of main beam bright-
ness temperature (Tmb)- This is related to via Tmb =
(Fes/Bcff) (cf. Downes 1989). The main beam efficiency, i?cff,
was 0.38 at 461 GHz and 0.5 at 330 and 352 GHz as ob-
tained by measurements of Saturn. The forward hemisphere ef-
ficiency. Fed, is 0.75 at 461 GHz and 0.9 at 330 and 352 GHz
(D. Muders, priv. comm.). The HHT is with an rms surface de-
viation of — 20/im (i.e. A/30 at 461 GHz) quite accurate. Thus
emission from the sidelobes should not be a problem.

Pointing was obtained toward Jupiter (continuum pointing)
and toward Orion-KL and R Cas (line pointing) with maximum
deviations of order 5". Observations were carried out in a posi-
tion switching mode with the off-position —1000" offset from
the source position.

I. Zinchenko et al.: HNCO in massive galactic dense cores

5

2.197 10^ 2/98 10^ 2.197 10^ 2.198 10^

Rest Frequency (MHz) Rest Frequency (MHz)

Fig. 2. Examples of 220 GHz HNCO spectra obtained at SEST
covering A' _ i = 0, 2 and 3 transitions

2.3. Data reduction and analysis

We have reduced the data and produced maps using the GAG
(Groupe d'Astmphysique de Grenoble) software package. The
measured spectra were fitted by one or more gaussian compo-
nents.

3. Results

3.1. One-point observations

HNCO was detected in 36 SEST sources (from 56 observed)
and in 22 OSO sources (from 27). Because of one source be-
longing to both samples, the total number of detected objects
is 57. In many cases K-i > transitions were detected too.
The gaussian line parameters are presented in Tables 5-14 (Ta-
bles 7-14 are available only electronically). It is worth noting
that a single-gaussian fit is clearly insufficient in many cases
because the lines have broad wings and other non-gaussian fea-
tures. Therefore, the values in the tables give only a rough rep-
resentation of the line profiles (the integrated intensities were
obtained by integrating over the lines in most cases).

Table i summarizes the 220 GHz SEST results for HNCO
10o,io - 9o,9 and C^^O. The Onsala 5o5 - 4o4 and C^^^O results
are presented in Table ^. The 220 GHz results for the A'_i — 2,
3 ladders ai-e given in Tables 7, 8. The 1 10 and 154 GHz SEST
data are displayed in Table 9. The Onsala 88 GHz data are sum-
marized in Table 10. The Effelsberg data are presented in Ta-
ble 11. Tables 12-14 contain the HHT data. We fitted the Ef-
felsberg spectra with 3-componentgaussians with fixed separa-
tions corresponding to the hyperfine structure of the Iqi — Oqo
transition.

Examples of measured spectra are given in Figs. ^ ^. Fig. ^
shows spectra of a few sources covering K^i = 0, 2 and 3
transitions at 220 GHz. Fig. ^ presents HNCO spectra in the
HNCO AT- 1=0 transitions at different wavelengths for sev-
eral sources.

The HNCO line profile in Orion KL can be decomposed
into at least two components which likely correspond to the so-
called classical "Hot Core" and "Plateau" outflow components
(see, e.g., Harris et al. 1995). The ratio between these compo-
nents is practically the same for the 5o5 — 4o4, 7o7 — 6o6 and
lOo.io — 9o,9 lines: ~ 60% of the emission originates from the
"Plateau" outflow source. The other lines do not allow such
decomposition due to their weakness or blending with other
spectral features.

An inspection of Table || shows that the derived C^'^O
velocities are systematically lower (more negative) than the
HNCO ones. The difference is — 1 km/s on the average. This
can be an instrumental effect: the C^^O line was located far
away from the center of the spectrometer band and a possible
non-linearity in the frequency response could lead to the appar-
ent displacement of the line on the velocity axis. This remark
is applicable also for the higher Ar_i HNCO lines.

3.2. Maps

In order to estimate source sizes and their spatial associa-
tion with YSO and infrared (IR) sources we mapped 2 south-
ern sources in the lOo.io — Qo.g HNCO line and Orion KL,
W49N and W51M in the 15o,i5 - 14o,i4, 16o,i6 - 15o,i5 and
21o,2i ^ 2Go,2o lines. W51M was mapped also in the 5o5 — 4o4
line. Three of these maps are presented in Fig. Q

The sources remain spatially unresolved. E.g. for
G 301.12-0.20 we obtain a FWHM « 29" in right as-
cension (from the strip scan across the map) which is very
close to the beam size at this frequency (24").

5.5. Detection of the K-i = 5 HNCO transition

The highest Ar_i HNCO transition reported so far was Ar_i =
4 (the IO4 - 94 Hne) in Orion (Sutton et al. 1985). This Une
is located on the shoulder of the strong C^^O J = 2 — 1 line.
In Fig. ^ we show parts of our Orion 220 GHz low resolution
spectrum and 461 GHz spectra with Ar_i = 2, 3, 4 and even 5
features (the K_i = 1 transition is outside our band). The rest
frequencies are assumed to be equal to those given in the JPL
catalogue for the strongest components of the corresponding
transitions (for A'_ 1 = 2 at 220 GHz we took the mean of the
frequencies of the two strongest components).

There is a weak bump in the redshifted C^^O J = 2 — 1
wing which can be attributed to HNCO IO4 — 94. Due to the
uncertainty in fitting the C^^O line profile the intensity of the
HNCO feature cannot be reliably determined but it is lower
than reported by Sutton et al. (1985). Our best estimate for the
integrated intensity is J Tj^hdv ^ 0.7 Kkm/s, but a reliable
error cannot be given.

There is also a feature at the A'_i = 5 frequency in the
220 GHz spectrum. It is located in the wing of a C2H3CN line.
The integrated intensity is J T^\^dv — 0.27 ± 0.06 Kkm/s.
The identification of this feature with HNCO seems to be reli-
able. The only other candidate is the C2H5OH 142i_2 — 13ii,2
line at 219391.81 MHz. However, there is no sign of other

6

I. Zinchenko et al.: HNCO in massive galactic dense cores

Table 5. C^^O (2 — 1) and HNCO (10o,io — 80,9) integrated intensities and gaussian line parameters at the indicated positions (cf.
Tables |l|,^ measured at SEST. The numbers in the brackets are the statistical uncertainties in the last digits (standard deviations).

C"0(2- 1) HNCO (lOo.io - 9o,9)

Source

Aa

/ Tmhdv

Vlsr

/ Tmbdv

Vlsr

A\/

(")

(")

(K-km/s)

(K)

(km/s)

(km/s)

(K-km/s)

(K)

(km/s)

(km/s)

G 261.64

20

17.22(05)

3.73(01)

13.76(01)

4.02(02)

0.63(05)

0.15(01)

10 z^ c\/ \ z^ \

13.69(16)

OT /

3.87(39)

G 264.28

-40

3.41(04)

f\ C\f\/ f\-\ \

0.90(01)

5.39(02)

3.36(05)

0.17(03)

0.07(01)

6.12(25)

-\ / ^ y' \

2.32(46)

G 265.14

-40

19.16(06)

4.63(02)

6.75(01)

3.68(02)

G 267.94

14.08(07)

2.59(01)

1.04(01)

4.68(03)

G 268.42

29.38(06)

5.57(01)

2.55(00)

4.60(01)

0.28(05)

0.07(01)

r\ 1 / ^ \

3.01(36)

r\ /T \

3.92(73)

G 269.16

40

24.95(08)

4.30(02)

c\ T"^ / r\r\\

9.72(00)

c f\c\ 1 r\'~\\

5.09(02)

A '7A/A'7\

0.79(07)

A 1 /A'^\

0.18(02)

A /' /I / 1 0\

9.64(18)

/I AC / d /I \

4.05(44)

'~\'~if\ '~\ /'

G 270.26

-20

20

16.28(05)

3.24(01)

8.67(01)

4.55(02)

1 AT/A/'\

1.07(06)

0.24(02)

A r^o / 1

9.92(12)

A 'O C /O A\

4.25(30)

G 285.26

10.16(06)

1.76(01)

-\ A*^ / fW \

2.47(01)

5.11(04)

0.50(09)

f\ f\r-i / f\-\ \

0.07(01)

1.64(65)

7.06(156)

G 286.20

40

-40

1 A C / r\ A \

14.53(04)

3.14(01)

-20.86(01)

A \ A / f\ 1 \

4.14(01)

'~\C\ 1 '~\^

G 291.27

^0

^0

28.39(06)

5.21(01)

-23.95(01)

4.93(01)

0.56(07)

0.12(02)

-23.80(25)

/I 'OC/'7A\

4.25(79)

G 291.57

7.59(07)

1.13(01)

13.12(03)

6.03(07)

G 294.97

15.77(05)

4.02(02)

-9.24(01)

001 /r\ I \

3.31(01)

G 300.97

40

23.13(05)

4.39(01)

-43.89(01)

A /' A / r\ I \

4.64(01)

G 301.12

80

-80

A A T^/10\

44.35(18)

7.13(03)

^O.Ol(Ol)

5.44(03)

4.15(08)

0.61(01)

T / r\z~ \

-39.37(06)

6.43(16)

G 305.20

22.98(07)

3.41(01)

^2.13(01)

6.39(02)

1 A / t A\

1.84(10)

AIT / A 1 \

0.17(01)

A f\ TO /'~\ /" \

-40.73(26)

1 A 1 'O /T 1 \

10.12(71)

G 305.36

28.45(08)

A A^ / fW \

4.47(01)

-36.60(01)

5.78(02)

Oi^o f\r\

G 308.00

13.38(05)

3.01(01)

-23.22(01)

3.85(02)

0.42(06)

0.14(02)

-21.63(16)

/I /o c \

2.42(35)

G 308.80

16.54(07)

2.73(01)

00 11 / r\ 1 \

-33.11(01)

5.27(03)

1 1 'O / AA\

1.12(09)

0.15(01)

-32.57(29)

7.13(62)

G 308.92

25.96(07)

5.07(02)

-51.46(01)

4.55(02)

A 1 / A A \

0.18(04)

A 1 A/AO \

0.10(03)

CI A0/10\

-51.08(18)

1.66(43)

G 309.92

22.08(07)

4.32(02)

f\/c\ 1 \

-57.60(01)

4.35(02)

oil' '1^

G 316.77

20

20

15.30(08)

2.51(01)

A 1 r\c /r\ 1 \

^1.05(01)

c /I /r\ A \

5.48(04)

G 316.81

20

17.51(08)

2.85(01)

-39.82(01)

5.63(03)

G 318.05

27.62(07)

5.63(02)

cf\ A /r\ 1 \

-50.48(01)

All /r\o\

4.31(02)

G 323.74

20

4.42(06)

1.11(02)

-DU.4o((JZ)

3.29(06)

A '~\r\

G ilA.lO

30

17.16(09)

2.39(02)

-89.27(02)

6.11(04)

G 326.47

12.13(09)

2.18(02)

-42.61(02)

4.59(05)

1.41(13)

A 1 /' /A'^\

0.16(02)

/II OA /O C\

-41.39(35)

1 T /r\c\\

8.17(99)

G 326.64

35.04(08)

7.97(02)

Ar\ A /f\f\\

-4U.z4(UU)

3.95(01)

A A £i{f\0\

0.46(08)

A \ r\(r\ \ \

0.10(01)

OA A A f\\

-39.24(40)

4.31(71)

G 328.24

24.60(09)

3.21(01)

An /' c\/r\ 1 \

-43.69(01)

7.11(03)

G 328.30

c 00/1 /^\

25.83(10)

3.56(02)

-92.92(01)

/' /IT / AO \

6.47(03)

G 328.81

55.86(15)

8.60(03)

-42.67(01)

5.34(02)

4.65(10)

A /'A/A^X

0.60(02)

At 0^7 /A'7\

^1.37(07)

Z" 0'7/^A\

6.37(20)

G 329.03

13.89(08)

2.11(01)

^4.19(02)

5.57(04)

1 T 1 / r\r\\

1.71(09)

0.28(02)

/lO Tr\/ic\

-43.70(15)

5.69(41)

G 330.88

43.65(23)

6.92(04)

-63.38(01)

5.45(04)

2.59(14)

0.34(02)

-62.78(18)

T /^ / ^ c \

7.08(45)

G 330.95

20

20

59.42(16)

6.79(02)

-92.16(01)

/ AO \

8.36(03)

O / 1

2.83(12)

A T7 /A 1 \

0.27(01)

A 1 00 /'^ 1 \

-91.38(21)

A T~l / A A\

9.77(44)

G 331.28

40

-20

21.07(09)

3.86(02)

-89.39(01)

5.01(03)

G 332.83

-20

85.39(13)

10.90(02)

-57.94(01)

7.13(01)

^ / 1 \

3.24(13)

r\ A A / /~\ \

0.44(02)

CT 00/1 /i\

-57.38(14)

y' /O \

6.92(33)

G 333.13

59.60(17)

8.28(03)

-52.55(01)

y' 1 /i~\o \

6.13(02)

1.07(13)

^ 1 'O / /~\ \

0.12(02)

-52.22(38)

6.89(140)

G 333.60

20

A f\ f\ A / -i A\

40.04(14)

4.10(02)

-49.60(02)

a ^ /i^ ^ "\

8.78(04)

G 337.40

20

20

43.56(10)

7.48(02)

-41.77(01)

^ //^ 1 \

5.29(01)

2.81(16)

^ T //^o \

0.37(03)

A r\ z^ c\/ t z^ \

-40.69(16)

G 340.06

34.59(13)

4.60(02)

-54.62(01)

6.59(03)

1.25(13)

f\ -l r-i / f\rw

0.17(02)

CO CC/OT\

-53.55(37)

6.97(81)

G 343.12

16.80(17)

2.84(03)

-31.08(03)

5.24(07)

0.59(11)

0.13(02)

-31.12(43)

4.34(87)

G 345.00

27.05(40)

2.74(06)

-27.58(07)

7.12(16)

3.56(13)

0.39(02)

-26.85(15)

c 1 /o \

8.51(38)

G 345.01

A A ^ / c\a\

44.61(08)

6.99(01)

-15.04(01)

5.79(01)

1.39(10)

0.25(02)

-14.27(18)

C /O / /I T \

5.28(47)

G 345.41

o /' r\c\/r\r\\

36.99(09)

6.46(02)

-22.11(01)

C IT /A'^\

5.17(02)

G 345.51

32.29(13)

5.03(03)

10 A C / C\ \ \

-18.45(01)

^ 1 //^ \

5.31(03)

1.10(12)

r\ 'O //~\o \

0.22(03)

-17.34(25)

4.70(60)

G 348.55

31.34(12)

4.32(02)

-16.16(01)

6.36(03)

G 348.73

76.64(12)

11.80(02)

-12.68(00)

5.78(01)

1.38(09)

0.25(02)

-11.16(17)

5.12(38)

G 350.10

27.58(11)

2.61(01)

-68.83(02)

9.98(04)

0.87(14)

0.07(01)

-68.18(99)

11.77(160)

G 351.41

57.97(20)

9.05(04)

-7.51(01)

5.57(03)

3.68(20)

0.42(03)

-7.60(18)

7.45(62)

G 351.58

24.04(12)

3.34(02)

-96.43(00)

6.37(04)

2.20(10)

0.36(02)

-94.72(13)

5.67(31)

G 351.78

88.86(38)

10.40(06)

-3.55(02)

6.11(04)

10.54(15)

1.08(02)

-2.95(05)

8.59(14)

G 353.41

55.52(10)

8.52(02)

-16.35(01)

6.05(01)

1.30(12)

0.25(02)

-16.11(22)

4.92(52)

G 359.97

20.16(09)

5.28(02)

17.71(01)

3.57(02)

Orion KL

53.52(135)

5.16(23)

7.69(11)

6.99(35)

52.50(70)

3.90(08)

7.16(07)

10.38(23)

SgrA

31.22(48)

1.50(03)

11.96(15)

17.51(37)

12.95(25)

1.07(02)

14.19(10)

10.96(26)

S231

8.10(09)

1.55(02)

-17.10(03)

4.71(07)

0.36(08)

0.15(03)

-15.45(26)

2.25(55)

S255

15.50(06)

3.50(02)

6.67(01)

3.97(02)

0.53(07)

0.15(02)

8.24(23)

3.26(47)

Fig. 3. HNCO K-i = lines in four selected sources. For the Iqi — Oqo transition a 3-component gaussian (according for
hyperfine structure) is superposed

8

I. Zinchenko et al.: HNCO in massive galactic dense cores

Table 6. C^^O (1 — 0) and HNCO (5o5 — 4o4) integrated intensities and gaussian line parameters at the indicated positions (cf.
Tables measured at Onsala. The numbers in the brackets are the statistical uncertainties in the last digits (standard deviations).

w l^i —

riLw^yj {DOS

— 404.)

Source

Aq

AS

/ Tmhdv

Vlsr

AV

Tmb

Vlsr

AV

(")

(")

(K-km/s)

(K)

(km/s)

(km/s)

(K-km/s)

(K)

(km/s)

(km/s)

G 121.30

5.94(18)

1.30(07)

-17.32(09)

3.55(22)

0.98(14)

0.40(06)

-17.57(16)

2.30(37)

G 123.07

3.19(14)

0.78(06)

-30.28(08)

3.81(32)

1.07(20)

0.12(02)

-31.14(84)

8.21(182)

G 133.69

-40

8.39(20)

1.22(03)

^2.38(08)

5.94(19)

0.55(14)

0.14(03)

-43.37(51)

3.77(106)

G 133.95

10.07(19)

1.71(04)

^7.60(05)

5.12(13)

1.38(24)

0.13(03)

-46.75(77)

9.95(231)

G 170.66

2.19(20)

0.56(05)

-15.27(24)

3.37(35)

G 173.17

4.50(27)

0.87(07)

-19.83(13)

4.78(40)

G 173.48

3.85(19)

0.76(04)

-16.22(12)

4.77(28)

1.25(20)

0.33(08)

-15.98(20)

2.68(69)

G 173.72

3.32(20)

0.97(09)

-16.83(12)

3.64(38)

0.49(14)

0.21(06)

-17.58(32)

2.22(72)

G 188.95

3.40(33)

0.70(09)

2.93(20)

4.54(66)

G 34.26

31.00(28)

4.45(04)

57.80(03)

6.47(07)

2.89(20)

0.54(04)

30.73(18)

5.55(41)

G 40.50

4.69(13)

0.81(03)

32.61(07)

5.22(20)

0.60(20)

0.14(08)

32.54(47)

3.94(237)

G 43.17

25.71(40)

1.45(02)

7.91(13)

14.84(28)

3.31(28)

0.24(02)

8.77(56)

12.97(118)

G 49.49

48.02(50)

2.99(04)

57.82(08)

12.45(19)

7.93(30)

0.69(03)

56.77(19)

10.20(50)

G 60.89

6.17(24)

1.33(07)

22.40(07)

4.35(26)

0.91(25)

0.14(03)

21.15(94)

7.07(143)

G 61.48

4.65(19)

1.12(06)

21.85(09)

3.82(23)

G 69.54

9.26(21)

1.98(05)

11.19(04)

4.52(14)

2.80(22)

0.43(04)

11.04(20)

5.38(57)

G 70.29

6.06(16)

0.70(02)

-24.28(10)

7.60(25)

G 75.78

8.96(18)

1.43(04)

-0.04(06)

5.28(14)

1.58(22)

0.16(02)

-1.88(66)

9.28(145)

G 77.47

5.59(16)

1.16(03)

0.50(07)

3.98(13)

1.00(20)

0.14(03)

0.98(72)

6.75(151)

G 81.72

19.07(22)

3.74(04)

-2.77(03)

4.49(06)

3.20(20)

0.53(04)

-2.71(19)

5.55(41)

G 81.77

8.23(13)

1.97(05)

-4.16(04)

3.65(11)

1.12(12)

0.25(04)

-3.97(19)

3.71(57)

G 81.87

13.45(24)

2.44(05)

9.36(05)

4.73(10)

2.37(28)

0.36(06)

9.52(31)

5.61(106)

G 92.67

4.22(16)

1.12(03)

-6.12(08)

3.66(12)

0.82(14)

0.21(04)

-5.26(37)

3.68(70)

G 99.98

5.21(14)

1.33(04)

0.67(06)

3.64(11)

0.46(10)

0.23(07)

0.67(15)

1.21(39)

G 108.76-0.95

6.79(16)

1.69(08)

-50.78(04)

3.77(20)

0.63(12)

0.23(06)

-50.63(22)

2.00(54)

G 108.76-0.98

11.97(17)

3.09(05)

-51.13(04)

3.51(06)

0.67(16)

0.16(04)

-51.08(45)

4.01(110)

G 111.53

12.22(22)

1.80(04)

-56.23(05)

6.18(14)

2.74(20)

0.55(05)

-56.44(16)

4.66(42)

100 80 60 20 10 -10 -20 -30 -40 20 10 -10 -20

Aa (") Aa (") Aa (")

Fig. 4. The HNCO lOoao - 9o,9 integrated intensity maps of G301. 12-0.20 and G270.26+0.83 and HNCO 21o,2i - 20o,20
integrated intensity map of Orion KL. The levels start from 15% of the peak intensities in steps of 7.5%. The peak intensities
equal to 4.4, 1.5 and 95.0 K kms^i for G301. 12-0.20, G270.26+0.83 and Orion KL, respectively. The beam widths are 24" for
the first two objects and 18" for Orion KL. For the first two objects the large crosses indicate IRAS positions, small stars show
NIR sources, triangles mark H2O masers, squares correspond to OH masers and diamonds show methanol masers (for references
see Lapinov et al. 1998). Open circles mark UC H II regions (Walsh et al. 1998). The IRAS uncertainty ellipses are shown. For
Orion KL only the IRc 2 position is indicated

I. Zinchenko et al.: HNCO in massive galactic dense cores

9

Velocity (km/s)

20
Velocity (km/s)

Fig. 5. Left panel: parts of the Orion low -resolution 220 GHz spectrum corresponding to higher K-i HNCO transitions. The
profiles are aligned in velocity. No baselines are subtracted but the subspectra are shifted along the j/-axis for clarity. The dashed
vertical line corresponds to Vlsr = 6.5 kms^^. Right panel: the same for the 461 GHz spectrum. Here the dashed vertical line
corresponds to Vlsr = 7.5 kms~^

ethanol lines in our spectrum so we reject this alternative. In
the 461 GHz spectrum the K^i — 5 feature is clearly detected.
Its integrated intensity is / T^\ydv ~ 1.4 ± 0.2 Kkm/s.

3.4. Hyperfine splitting, HN CO and optical depths

The HNCO lines are split into several hyperfine components
mainly due to the ^"^N spin. This splitting is clearly seen in the
loi - Ooo transition (Fig. |) at 22 GHz. EarUer HNCO hyper-
fine structure in the Iqi — Oqo line was only observed in the dark
cloud TMC-1 ( Brown 198 1[ ) where possible deviations from
the optically thin LTE (Local Thermodynamic Equilibrium) in-
tensity ratios (3:5:1) were found. In our spectra the hyperfine
ratios are consistent with the optically thin LTE values. Taking
into account the measurement uncertainties, an upper limit on
the optical depth in this transition for the sources detected in
Effelsberg is t ^ 1 .

To the best of our knowledge no isotopomer of HNCO ex-
cept the main one has been detected in space yet. This de-
tection would be important for estimates of HNCO optical
depths which are believed to be small (e.g. Jackson et al. 1984,
Churchwell et al. 1986). The frequency separations between

the HN^'^CO and the main isotopomer lines are rather small
corresponding to a few km/s, so in sources with broad lines like
Orion A or Sgr A the HN^^CO fines will be blended. However,
there are some strong HNCO sources in our sample with nar-
rower lines which show features attributable to HN^'^CO. The
most reliable one is seen in the G 301.12-0.20 10o,io ~ 9o,9
spectrum (Fig. ||). A weak feature on the blue shoulder of the
main isotope line is very close in frequency to the expected
location of the HN^^CO line .

For comparison we show in addition to HNCO also the
C'^'^S spectrum. It is noteworthy that there is no bump in this
spectrum corresponding to the discussed feature in HNCO.

The fine we identify with HN^^CO is shifted by 0.65 ±
0.21 MHz from the expected HN^^CO transition frequency.
This 3(7 shift, if it is significant, cannot be explained by instru-
mental effects like in the case of our C^^O data because the fea-
ture is very close to the main isotope line. The shift greatly ex-
ceeds the uncertainty of the transition frequency derived from
the laboratory data (Winnewisser et al. 1978) which is 25 kHz.
This makes the identification questionable. Detection of other
HN^'^CO lines would be important in this respect. There is no
corresponding feature in the 7o7 — 6o6 HNCO spectrum (the

10

I. Zinchenko et al.: HNCO in massive galactic dense cores

Q.

E

0)

0.5 -

E

o

X)

-60 -40

Velocity (km/s)

-20

Fig. 6. The HNCO 10o,io - 9o,9 line in G301. 12-0.20 in com-
parison with the C'^'*S(2-1) line. The expected location of the
HN^^co 10o,io - 9o,9 line is shown

5o5 — 4o4 spectrum is too noisy). This could mean that the op-
tical depth in this transition is significantly lower Indeed, at
sufficiently high temperatures (> 30 K) it can be about 2 times
lower than in the 10o,io — 9o,9 transition according to Eq. O)

(see the discussion in Sect. 4.2)

If our identification of the discussed line with HN^'^CO is
correct we can estimate the optical depth assuming the same
excitation as for the main isotopomer For G 301.12-0.20 we
obtain r(HNCO) « 15 if we assume the terrestrial ^^C/^^C
isotope ratio C^C/^^C = 89) and « 7 for ^^c/i^C = 40. A high
optical depth in the 10o,io— 9o,9 HNCO line does not contradict
our conclusion of low optical depth in the Iqi — Ooo transition
because the line strengths for these transitions are different (see
discussion in Sect. 4.2). Therefore, the optical depth in some
lines of the main isotopomer might be rather high. This contra-
dicts the usual assumption of low optical depth in all HNCO
lines (e.g., Jackson et al. 1984, Churchwell et al. 1986) and
could imply serious consequences for the analysis of HNCO
excitation and abundances.

4. Discussion

4.1. Comparison with other HNCO data

Most of our HNCO sources are new detections. Only few were
included in the surveys of Jackson et al. (1984) and Church-
well et al. (1986). A direct comparison with the intensities
measured by Jackson et al. is impossible due to different tem-

perature scales. Common detected sources are Orion KL and
W51. Their upper limit for W3(OH) does not contradict our
value if we take into account the difference in the temperature
scales. The upper limits for the Iqi — Ooo transition obtained
by Churchwell et al. do not contradict our results taking into
account the differences in the beam sizes and efficiencies.

As mentioned above, towards Orion KL several HNCO
lines were observed at 220 GHz by Sutton et al. (1985). Their
results agree in general with our measurements though there is
a discrepa ncy concerning the intensity of the K^i — 4 transi-
tion (Sect. |33|).

It is worth noting that while at 22 GHz and at 110 GHz
(as obtained by Jackson et al. 1984) the brightest source of
HNCO emission is the Galactic center, at 220 GHz the situation
changes and Orion becomes the brightest source with several
other sources approaching Sgr A in intensity. Apparently this
is caused by differences in excitation.

4.2. Rotational diagrams

As a first step in the excitation analysis we construct traditional
rotational diagrams for our sources. For a recent discussion of
this method see e.g. Goldsmith & Langer (1999). This means
a plot of the column density (N^) per statistical weight (gu)
of a number of molecular energy levels, as a function of their
energy above the ground state (i?u)- In local thermodynamic
equilibrium (LTE), this will just be a Boltzmann distribution,
so a plot of ln{Nu/gu) versus Eu/k will yield a straight line
with a slope of 1/Tr. The temperature inferred is often called
the "rotational temperature".

Actually from the measurements we do not obtain directly
the column densities. The measured quantity is the line inten-
sity. In an optically thin case for Tex >> Jbg (Tox is the exci-
tation temperature of the transition and Tbg is the background
temperature)

log

= log

N

kTT,

lege

(1)

where W is the integrated line intensity, /b is the beam dilution
factor, S is the line strength, ji^ is the appropriate component
of the dipole moment, N is the total column density and Q{T'q)
is the partition function.

The quantity on the left hand side of Eq. (|l|) can be derived
from the molecular data. Plotting it versus E^ we can find the
rotational temperature (from the slope) and the total column
density (from the intercept).

Some problems can arise from an uncertainty in the beam
filling factor As shown in Fig. Qthe sources are probably unre-
solved. Assuming that the source size is the same for all HNCO
transitions in a given source and that the source size is small
with respect to the beam, we reduced all data to the same beam
size, the SEST HPBW at 220 GHz, i.e. 24".

For Orion the highest observed transition lies ~ 1300 K
above the ground level. For other sources we managed to ob-
serve transitions up to ~ 450 K above the ground state. Exam-
ples of the rotational diagrams are presented in Figs. H.

I. Zinchenko et al.: HNCO in massive galactic dense cores

11

CO

\ 100 ?00 100 700 300 400

20 40 60 100

Eu (K)

Fig. 7. Rotational diagrams for selected sample sources (W =
/ Tynhdv, S is the line strength). Filled squares correspond to
the measured values and the open squares to the values cor-
rected taking into account the beam sizes (see text). The dia-
mond on the Sgr A plot corresponds to the data from Lindqvist
et al. 1995

The measured integrated intensities are represented by
filled squares (/b — 1). The corrected results are plotted by
open squares m Figs. 0, H One can see that they much better
correspond to each other than the uncorrected values.

The rotational diagram for Orion is presented in Fig. |^. The
rotational temperature from this plot is Tiot ~ 25 K for the
lowest transitions and Tiot ~ 530 K for the highest transitions.
The latter one is a very high value even for Orion KL. But in
principle the diagram shows a range of rotational temperatures.
We represent it by 3 components as shown in Table [l^. A sepa-
rate fit to the K-i — transitions gives T^ot « 80 K (although
this fit is not very satisfactory).

The rotational temperatures and column densities derived
from rotational diagrams are summarized in Table |l5|. In this
analysis we assume that the sources are optically thin in the
observed transitions. This contradicts the tentative detection of
HNi^CO in G 301.12-0.20. The effects of high optical depth
on rotational diagrams have been analyzed recently by Gold-
smith & Langer (1999). In optically thick case the column den-
sity in the upper level of the transition (iVu) is underestimated
by the factor of t/(1 — e~'^) and, therefore, corresponding
points in the population diagram lie lower than they should. In
general, for linear molecules it produces a curvature resembling
that seen in the diagrams for Orion and some other sources. It
is caused by the fact that the optical depth exhibits a peak for
transitions with the excitation energy ^ kT (Goldsmith
& Langer 1999). However, for nonlinear molecules the opti-
cal depth effect rather leads to a "scatter" in the population di-
agram, because transitions with significantly different optical
depth can have similar excitation energies.

500 1000

Eu (K)

Fig. 8. Logarithm of integrated line intensity divided by the line
strength and frequency versus upper-state excitation energy for
Orion A. The filled squares correspond to the measured values
and the open squares represent the values corrected for beam
width ratios (see text). The points corresponding to K^i =
transitions are encircled

There is a strong argument against high optical depth at
least for transitions with ~ 200 — 400 K in Orion. In
this range transitions with similar energies of the upper state
but with very different frequencies (belonging to different K-i
ladders) were observed. It is easy to estimate the expected ratio
of peak optical depths in the lines which is

For \E^ — E'^l << kT the exponential factor is close to unity.

In our data there are pairs of transitions with similar up-
per state energies. The 21o,2i — 20o,20 and 102,9/8 — 92,8/7
transitions have similar E^ ^ 200 K. However, the first one
has higher line strength and higher transition frequency; there-
fore, according to Eq. (||) it should have higher optical depth
than the second one. Then, it should be stronger influenced
by possible optical depth effects and the corresponding point
in Fig. H should lie lower than the point corresponding to the
102,9/8 ~ 92,8/7 transitions. However, this is not a case. Ac-
tually, the points are very close to each other and perhaps
slightly shifted in the opposite sense. The same is true for
the 212,20/19 ~ 202,19/18 and 103,8/7 ~ 93,7/6 transitions with
Eu ^ 400 K. We conclude that the optical depth for Orion in
these transitions should be low. Perhaps in some other transi-
tions or in other sources optical depths are as high as indicated
by our tentative HN^^CO detection. There is however no reason
to apply optical depth corrections to the bulk of our sources.

Transitions with low E^/k values are fitted by rather low
temperature models, ^ 10 — 30 K. Transitions between
higher excited states are related to higher rotational tempera-
tures up to Trot ^ 500 K. In Table |T5| we also present estimates
of the HNCO relative abundances. The hydrogen column den-
sities have been calculated from the C^^O data under the as-

12

1. Zinchenko et al.: HNCO in massive galactic dense cores

Table 15. HNCO rotational temperatures, column densities and
relative abundances.

rot

vein ;

Orion A

25

14.87

-8.06

150

15.00

530

14.61

G 301.12

24

13.73

-9.16

76

13.94

G 305.20

102

13.69

-9.33

G 308.80

236

13.89

-8.99

G 329.03

60

13.89

-8.91

G 330.88

133

13.96

-9.34

G 332.83

98

13.94

-9.65

G 337.40

75

13.72

-9.58

G 340.06

130

13.64

-9.06

550

13.98

G 345.00

88

13.93

-9.08

G 351.41

93

13.93

-9.17

320

13.96

G 351.58

70

13.65

-9.39

G 351.78

155

14.64

-8.88

SgrA

12

14.95

-8.20

S 158

28

14.03

-9.23

S255

170

13.68

-9.17

W49N

100

15.00

-8.58

W51 M

9

14.44

-8.74

38

14.84

73

14.49

W75 N

46

14.38

-8.92

W 75(OH)

37

14.38

-9.07

sumptions of LTE and a C^^O relative abundance of 1.7 10^^
(Frerking et al. 1982). Typical HNCO abundances are ~ 10~^.
Sgr A does not look very exceptional here. The relative HNCO
abundance in Sgr A is about the same as in Orion but the ro-
tational temperature is much lower. In contrast to many other
sources there is no high excitation temperature component in
Sgr A, indicating that the dense gas is probably cool. This
agrees with results from Hiittemeister et al. (1998) based on
SiO and C^^O. The opposite scenario, a hot highly subther-
mally excited low density gas component (n(H2) ^ 10^ cm^^^)
as observed by Hiittemeister et al. (1993) in ammonia toward
Sgr B2 is less likely, due to the correlations between HNCO
and SiO that will be outlined in Sects, i.4 and 4.6.

It is important to emphasize that our estimates give
lower limits to the relative abundance X(HNCO) =
7V(HNCO)/iV(H2) for at least two reasons. First, the HNCO
sources are much more compact than their C^*0 counterparts
and tend to be spatially unresolved. Our estimates give beam
averaged values and "real" abundances in regions of HNCO
line formation should be significantly higher Second, if the
HNCO optical depth is high we would underestimate its
column densities.

-8 F'

-8.5 -

-9 -

-9.5 -

AV(HNCO) (km/s)

Fig. 9. HNCO relative abundance versus the HNCO line width
for sources observed at SEST (filled squares) and in Onsala
(open squares)

Next, we have to mention that all these estimates refer to
the bulk of the cores. In the high velocity gas the HNCO abun-
dances are apparently much higher

One might think that better estimates of HNCO abun-
dances can be obtained from comparison with the dust emis-
sion rather than with C^'^O. As shown, HNCO probably arises
in "warm" environments and in the dust emission we see pref-
erentially a high temperature medium while in C^^O the re-
verse is true. However, interferometric observations in Orion
dBlake et al. 19961) show that HNCO K^i = 2 and dust distri-

butions do not entirely coincide. At the same time, as shown in

Sect. 4.5 , there is a tight correlation between the FIR emission
at 100 fim and C^®0(2-1) integrated line intensity. Therefore,
no large differences between estimates of HNCO abundances
by both methods can be expected. There are detailed studies
of dust emission towards some of our sources with comparable
angular resolution. E.g. Henning et al. (2000) show that total
gas column densities derived from dust and from C^^O(2-l) in
G301.12-0.20 coincide within a factor of 3.

In Fig. ^ we plot the HNCO abundances versus the HNCO
line widths. There is a trend of increasing the HNCO abun-
dance with increasing HNCO line width. This shows that the
HNCO production can be related to dynamical activity in the
sources.

Table |l5| and Fig. ^ indicate that abundances derived for the
sources which belong to the inner and to the outer Galaxy, re-
spectively, are about the same. Therefore, there is no significant
galactic gradient in HNCO abundance.

4.3. Physical conditions in regions of HNCO emission

Now we shall try to understand the physical conditions in re-
gions of HNCO emission detected by us. An important ques-
tion to start with is which excitation mechanism dominates, ra-
diative or collisional? And which gas parameters are implied
by each of them? To answer these questions properly would

I. Zinchenko et al.: HNCO in massive galactic dense cores

13

require a numerical model taking both into account. Useful
conclusions can, however, also be obtained by semi-qualitative
consideration presented below. We concentrate here on Orion
KL as the best studied source.

At first, we need an estimate for the size of the HNCO emis-
sion region. Our map presented in Fig. Qgives an upper limit of
^ 10" for the 21q 21 — 2O0.20 transition. Interferometric results
dBlake et al. 199^ give a size of ~ 2" x 4" for the X_i = 2
transition at 220 GHz. This can be probably considered as an
upper limit also for higher K^i ladders. On the other hand we
can obtain a lower limit on the source size from the comparison
of the brightness and excitation temperatures. For Tex ^ 500 K
(as follows from the population diagram) we obtain that the
lower limit on the beam filling factor for the K^i — 5 transi-
tions in Orion is ^ 2 lO^'^. Therefore, the effective size of the
emitting region is ^ 1" or 0.002 pc, i.e. 7 10^^ cm.

Let us consider the physical requirements in the case of
collisional excitation. The critical densities defined as rtc =
Au\/Cu\ (Aui is the spontaneous decay rate and Cui is the col-
lisional de-excitation rate; Scoville et al. 1980) are
for the IOq.io — 9o.9 transition and ^ 10^ cm^^^ :
20o,20 transition ( the collisional rates are ^ 10^^" s^^cm^-^ as
o btained from Sheldon Green's program available on Internet

10*^ cm"

http://www.giss.nasa.gov/data/mcrates4 ). Much higher densi

ties are needed for excitation of the transitions in the K^i >
ladders. This is caused by fast 6-type transitions between dif-
ferent K^i ladders. E.g. the spontaneous emission rate from
the K-i — 5 ladder to the K^i = 4 ladder is ~ 5 s~^. This
implies a critical density of ~ 10^^ cm^''. The gas kinetic tem-
perature should be <; 500 K.

Such conditions cannot be excluded. Walker et al. (1994)
derived from observations of vibrationally excited CS n <;
10" - 10^2 cm-3 and T ^ 1000 K in a region - lO^^ cm
from the stellar core toward IRAS 16293-2422. The question
is whether the required amount of such gas is consistent with
the observations.

Taking into account the lower limit on the source size the
mass of the hot dense gas (n ~ 10^^ cm^'^, T ^ 500 K) would
be > 100 Mq. Estimates of the hot core mass from dust con-
tinuum measurements give values of ^ 5 — 40 Mq (Masson
& Mundy 1988, Wright et al. 1992). Taking into account the
uncertainties in our estimations we cannot entirely exclude the
possibility of collisional excitation even for the K^i = 5 lad-
der but this appears to be an unlikely scenario.

For the lower K^i ladders the density requirements can be
significantly relaxed. E.g. for the 6-type transitions from the
K^i = 3 to the K^i = 2 ladder the spontaneous decay rate is
^ 1 s^^ and the critical density is ~ 10^*^ cm^'^.

The transitions in the K^i — ladder, of course, will be
also excited in this hot dense gas. However, the emission in
these lines will be dominated by a more extended lower density
component.

Now let us turn to radiative excitation. It requires sufficient
photons at the wavelengths corresponding to the 6-type transi-
tions between different K^i ladders, from ~ 300 to ~ 30 /im.
If the dilution factor is close to unity we need an optical depth

r <; 1 and a radiation temperature Tr ^ 500 K at least at
30 fim. As an upper limit to the source size we can take the
mean interferometric value of ~ 3". However, what will be
the IR flux and luminosity of such a source? For the flux at
30 Aim we obtain F(30;um) « 3 10'*(6's/l")2 Jy. The observa-
tional value is ^ 5 10^ Jy (van Dishoeck et al. 1998). There-
fore, the angular source size should be 9^ ^ 1'.'5 and the linear
size <} 10^^ cm. This practically coincides with the lower limit
on the source size derived from the beam dilution (see above).
Taking the dust absorption coefficient of k„i ^10^ cm^/g (Os-
senkopf & Henning 1994) we conclude that the gas density in
this region should be n 3 10^ cm^'^. In this case we have
no problem to reconcile the mass estimates with the available
data.

However, at longer wavelengths the IR pumping from such
a source might be not sufficient. Say, for r oc the optical
depth at 300 fim will be only 0.01. Therefore, we need even
higher gas column and volume densities and/or larger source
sizes at longer wavelengths. The latter implies the presence of
a temperature gradient in the source which is natural for an
internally heated object. The lower K^i ladders are apparently
excited by radiation with a lower effective temperature.

To conclude, it is much easier to explain the excitation of
the higher K^i ladders by the radiative process. The source
size in Orion should be 1" — 2" which agrees with the interfer-
ometric image in the K^i — 2 transition at 1.3 mm (Blake et
al. 1996).

The emission in the K^i = ladder should be more ex-
tended. For Orion again from a comparison between the bright-
ness and excitation temperatures the source size should be
^ 6". Such a large source size for the A'_i = transitions
implies that the radiative excitation via K^i > ladders will
become inefficient. Therefore, for the K^i = ladder col-
lisional excitation may dominate which implies gas densities
n ^ 10^ — 10' cm^'^. This scenario is supported by several
sources where the HNCO emission peak is significantly dis-
placed from any known IR source. The most obvious example
is G 270.26h-0.83 (Fig. Q). This implies either the presence of
a very dense prestellar core or a highly obscured young stellar
object at this location.

4.4. Comparison with C^^ O, CS and SiO data

An obvious step ahead to understand the properties of interstel-
lar HNCO emission is to compare our results with data from
other better studied species. The most reliable comparison can
be done with our C^**0 data which were observed simultane-
ously with HNCO.

Fig. [1^ shows a noticeable correlation between the HNCO
and C^'^O integrated line intensities. However, it is produced
apparently by the correlation between the line widths since the
correlation between HNCO and C^'^O peak line temperatures
is rather weak.

The plot of AT/(HNCO) versus AV{C^^O) looks rather
interesting. Concerning the 220 GHz transitions for the nar-
rowest C^^O lines the HNCO line width is smaller than that

14

I. Zinchenko et al.: HNCO in massive galactic dense cores

; "IE

; ' ■ ""1 ^:

; a I

- ' , ^

- ■ -

■ ■ .

■:

\' "

- ■

10 10
(C '^O) (K krr/s)

10 10
i(SiO) (K km/s)

o
o

10

10

o
o

b '
_ ■ ■ _

mm-

■ ■ , * . .

^ ■ ■

- ■ —

;^ 1 ■ , ■ , 1 E

- ■'-

_ ■ _
□ ■

□ ■ , !■

v."

- ■ -

■

1 1 1

T„b(c''0) (K) WSiO) (K)

- I ' -

■ ^

_ ^ -' ■
' ■ ■ ■

- , , 1 , 1 , , , 1 , , 1 , , -

- ' ' ' l; ' ' ' 1 -

/

- p' ■
/ . ■

/ □ a

r : ;

4 6 8 10
AV(C^*^0) (km/s)

10 20 30

AV(SiO) (km/s)

Fig. 10. HNCO IOq.io ~ 9o,9 integrated line intensities, peak
main beam temperatures and line widths versus correspond-
ing C^^O J — 2 — 1 and SiO J — 3 — 2 peak temperatures
for the SEST sample. Open squares correspond to those SiO
data which were obtained at slightly different positions than
HNCO. The dashed lines in the panel (c) correspond to equal
line widths of the compared species

of C^^O. With increasing C^^O linewidth, however, the HNCO
lines broaden faster and become broader than the C^^O lines.
An exception is Sgr A (not shown in the plot) but its C^^O
spectrum is strongly distorted by emission from the reference
position.

A similar comparison with the CS(2-1) data from
Zinchenko et al. (1995, 1998) and Juvela (1996) (not shown
here) shows even lower correlations between the line parame-
ters than in the case of C^*0. However, in this case the beam
sizes for CS and HNCO are different and even the central po-
sitions not always coincide.

In contrast, much better correlations exist between the
HNCO and SiO line parameters (the latter ones are taken from
Harju et al. 1998). Good correlations exist for both integrated
and peak intensities. The correlation between the line widths is
somewhat worse but one should take into account that the SiO
line widths were derived from the second moments of the line
profiles while the HNCO widths represent results of the gaus-
sian fits. Anyway, the correlation does exist and the SiO lines
are almost always broader than the HNCO lines.

Velocity (km/s)

Fig. 11. The HNCO IOq.io - 9o,9 Hne in Orion (thick line)
in comparison with the C^®0(2-1) (dotted Hne) and SiO(2-l)
(thin solid line). The latter ones are scaled to the same peak
intensity as HNCO 10o,io - 9o,9

A more detailed comparison with other species should in-
clude the line profiles. For Orion, such a comparison is dis-
played in Fig. ^ It shows that HNCO lines possess an extra
wing emission which is less pronounced than in SiO. A similar
picture is seen in some other sources.

This comparison shows that HNCO is closely related to
SiO which is thought to be produced primarily in shocks and
other energetic processes. The comparison with the presum-
ably optically thin 0^*^0(2-1) Une shows that the HNCO/CO
abundance ratio is apparently enhanced in high velocity gas
although to a lower degree than for SiO. Since the CO abun-
dance is usually assum e d to be constant in bipolar flows (e.g.,
ICabrit & Bertout 1992| ; [Shepherd & Churchwell E. 1996| ) we
see that HNCO abundances are enhanced relative to hydrogen,
too.

It is interesting to note that the interferometric data for
Orion (Blake et al. 1996) show that the spatial distributions of
SiO and HNCO are rather different. However, this does not ex-
clude a common production mechanism. E.g. these species can
be formed at different stages in the postshock gas.

4.5. Comparison with IR data

The correlation between HNCO integrated line intensities and
FIR flux, e.g. at 100 fim taken from IRAS data (Fig. looks
rather similar to the relationship between HNCO and C^^O
(Fig. |l^. This is natural because there is a rather tight cor-
relation between the 100 /im flux and the C^^O integrated line
intensity (Fig. [T^. Such a good correlation shows that C^*0
relative abundances are rather constant and justifies the usage
of the HNCO/C^^O ratio for estimation of HNCO abundances.

I. Zinchenko et al.: HNCO in massive galactic dense cores

15

10 10
F(100/im) (Jy)

Fig. 12. HNCO lOo.io — 9o,9 integrated line intensities versus
the FIR flux at 100 /im for the SEST sample. The open squares
correspond to the cases where there is a large (Jl HPBW) dis-
placement between the position observed in HNCO and the
IRAS position

10

10 10
F(100/im) (Jy)

Fig. 13. C^^O(2-l) integrated line intensities versus the FIR
flux at 100 /im for the SEST sample. The open squares corre-
spond to the cases where there is a large HPBW) displace-
ment between the position observed in C^^O and the IRAS po-
sition

4.6. HNCO chemistry

In the early work of Iglesias (1977) HNCO was suggested
to form via ion-molecule reactions. The sequence leading to
HNCO via electron recombination of H2NCO+ is initiated by
the formation of NCO+ (either by a reaction between CN and
O2 or between He+ and NCO; see also ^rown 1981 ). The pre-
dicted HNCO abundances from this reaction scheme are low.
The steady state fractional abundance is of the order 10^^° for
a model with n^i^ ^ lO** cm^^^ (Iglesias 1977), and still lower

The abundances derived from ion-molecule chemistry are
in contradiction with the observations, especially when HNCO
is believed to trace high density gas. Recently, a new neutral
gas-phase pathway has been suggested by Turner et al. (1999)
for translucent clouds: CN -I- O2 NCO + O followed by
NCO + H2 ^ HNCO + H. The importance of these reac-
tions can however be questioned, since 1) the abundance of
O2 in the interstellar space is poorly known; and 2) the sec-
on d reaction probably has an activation barrier of about 1000
K ( [Turner et al. 1999^ .

Chemistry models predict high fractional O2 abundances
(up to 10^^) at late stages of chemical evolution in

dense cores and in postshock gas (e.g Caselli et al. 1993

Bergin et al. 1998). However, the upper limits derived from ob-

servations towards several GMC cores (most recently by the
SWAS satellite; [Melnick et al."T999[ ) are about 10 ^, which in-
dicates that the oxygen chemistry is not well understood, yet.
O2 is destroyed by UV radiation and in po werful shocks (wit h
shock velocities greater than 26 kms^^; Bergin et al. 1998 ),

for higher densities, because the fractional ion abundances are
roughly inversely proportional to the square root of the gas den-
sity.

and is therefore likely thriving in relatively quiescent dense gas
or in regions associated with low velocity shocks. The same
should be true for HNCO if the reaction suggested by Turner et
al. (1999) is relevant.

The observed correlation between SiO and HNCO inte-
grated line intensities indicates the prevalence of shocks in the
HNCO emission regions. Shock heating can therefore provide
the means of overcoming the energy barrier in the reaction be-
tween NCO and H2, and thereby intensify the HNCO produc-
tion. On the other hand, the fact that the HNCO line widths are
smaller than those of SiO could be understood by the destruc-
tion of O2 in high velocity shocks.

In the light of the present observations the neutral reactions
suggested by Turner et al. (1999) appear to provide a plausible
production pathway of HNCO also in warm GMC cores. The
formation of HNCO via grain surface reactions, e.g. through
the desorption and subsequent fragmentation of some more
complex molecule is an alternative, which to our knowledge
has not yet been investigated.

5. Conclusions

We have presented the results of an HNCO survey of high mass
star-forming cores at frequencies from 22 to 461 GHz. The
main conclusions are the following:

1. HNCO is widespread in dense cores forming high mass
stars. The detection rate was ^ 70%. There is no signifi-
cant galactic gradient in its abundance as indicated by the
fact that abundances derived for the sources which belong
to the inner and to the outer Galaxy, respectively, are about
the same.

2. Transitions in higher K^i ladders, up to K^i = 5, are
detected. The excitation energy reaches ~ 1300 K above
the ground level.

3. HN^^CO is tentatively detected towards G 301. 12-0.20.
This implies an optical depth in the HNCO 10o,io — 9o,9

16

I. Zinchenko et al.: HNCO in massive galactic dense cores

line ^ 10 in this source. The optical depth in the Iqi — Oqo
transition is r ^ 1 for the sources detected in this line as
inferred from the hyperfine ratios.

4. The sources are compact with sizes ^ 20".

5. HNCO rotational temperatures vary from ~ 10 K to ^
500 K. Typical relative abundances are ~ 10~^. These in-
crease with increasing velocity dispersion.

6. The emission in the K-i > ladders is best explained by
FIR radiative excitation. In order to provide a sufficiently
large dust optical depth at FIR wavelengths taking into ac-
count the limitations on the source size, the gas density
should be n ^ 3 lO'^ cm"^; a temperature T ^ 500 K
is needed to excite the K^i = 5 emission in Orion KL.
The K-i = transitions can be colUsionally excited. The
required densities are n ^ 10^ — 10^ cm~^.

7. HNCO correlates well with SiO and does not correlate with
CS which is a typical high density probe. HNCO abun-
dances are enhanced in high velocity gas. Probably HNCO
production is related to shocks as for SiO. A plausible path-
way is gas-phase neutral-neutral reactions at high (> 1000
K) temperatures to overcome an activation barrier that is
likely inhibiting the NCO + H2 — > HNCO -|- H reaction in
a cool interstellar medium.

Acknowledgements. We are very grateful to Dr. J. Harju for his con-
tribution to this work, to Dr. Lars E.B. Johansson for the help with the
observations in Onsala, to the SEST staff, to Alexander Lapinov for
calculating HNCO line strengths and to the referee. Dr. CM. Walm-
sley, for the very useful detailed comments. l.Z. thanks the Helsinki
University Observatory and Max-Planck-Institut fiir Radioastronomie
for the hospitality. He was also supported in part by the DFG grant
436 RUS 1 13/203/0, INTAS grant 93-2168-ext, NASA grant provided
via CRDF RPO-841 and grants 96-02-16472, 99-02-16556 from the
Russian Foundation for Basic Research. This research has made use
of the Simbad database, operated at CDS, Strasbourg, France.

References

Baars, J., Martin, R.N., 1996, Rev. Mod. Astron. 9, 111
Bergin E. A., Melnick G.J., Neufeld D.A. 1998, ApJ 499, 777
Blake G.A., Mundy L.G., Carlstrom J.E., et al., 1996, ApJ 472, L49
Brown R.L. 1981, ApJ 248, LI 19
Cabrit S., Bertout C, 1992, A&A 261, 274
Caselli R, Hasegawa T.I., Herbst E. 1993, ApJ 408, 548
Churchwell E., Wood D., Myers RC, Myers R.V. 1986, ApJ 305, 405
Dahmen G., Huttemeister S., Wilson T.L., Mauersberger R., 1997,
A&AS 126, 197

Downes, D., 1989, Evolution of Galaxies - Astronomical Obser-
vations, Lecture Notes in Physics 333, eds. 1. Appenzeller, H.
Habing, P. Lena, Springer Verlag, Berlin, p353
Frerking M.A., Langer W.D., Wilson R.W., 1982, ApJ 262, 590
Goldsmith RF., Langer W.D., 1999, ApJ 517, 209
Groesbeck T.D., Phillips T.G., Blake G.A., 1994, ApJS 94, 147
Harju J., Lehtinen K., Booth R., Zinchenko 1., 1998, A&AS 132, 21 1
Harris A.I., Avery L.W., Schuster K.-F, Tacconi L.J., Genzel R., 1995,
ApJ 446, L85

Henning Th., Lapinov A., Schreyer K., Stecklum B., Zinchenko I.,

2000, A&A, submitted
Hiittemeister S., Wilson T. L., Henkel C, Mauersberger R., 1993,

A&A 276, 445

Hiittemeister S., Dahmen G., Mauersberger R., et al., 1998, A&A 334,

646

Iglesias E. 1977, ApJ 218, 697

Jackson J.M., Armstrong J.T., Barrett A.H. 1984, ApJ 280, 608

JuvelaM., 1996, A&AS 118, 191

Kuan Y.-J., Snyder L.E. 1996, ApJ 470, 981

Lapinov A.V., Schilke P, Juvela M., Zinchenko 1., 1998, A&A 336,
1007

Lindqvist M., Sandqvist Aa., Winnberg A., Johansson L.E.B., N3Tnan

L.-A. 1995, A&AS 113, 257
Martin-Pintado J., de Vicente P., Wilson T.L., Johnston K.J., 1990,

A&A 236, 193
Masson C.R., Mundy L.G., 1988, ApJ 324, 538
Mehiick G.J., Stauffer J.R., Ashby M.L.N, et al. 1999, BAAS 194,

4708

Ossenkopf V., Henning Th., 1994, A&A 291, 943
Ott M., Witzel A., Quirrenbach A., et al., 1994, A&A 284, 331
Sato R, Hasegawa T., Whiteoak J.B., Shimizu M., 1997, In: lAV Sym-
posium No. 184. "The Central Regions of the Galaxy and Galax-
ies". Kyoto, Japan, 17-30 August, 1997, p. 98
Schulz, A., Henkel, C, Beckmann, U., et al., 1995, A&A 218, 24
Schilke R, Groesbeck T.D., Blake G.A., Phillips T.G., 1997, ApJS
108, 301

Scoville N.Z., Krotkov R., Wang D., 1980, ApJ 240, 929
Shepherd D.S., Churchwell E., 1996, ApJ 472, 225
Snyder L.E., Buhl D., 1972, ApJ 177, 619

Sutton E.G., Blake G.A., Masson C.R., Phillips T.G., 1985, ApJS 58,
341

Townes C.H., Schawlow A.L. 1975, Microwave Spectroscopy, Dover

Publications, New York
Turner B.E., Terzieva R., Herbst E. 1999, ApJ 518, 699
van Dishoeck E.F., Wright CM., Cemicharo J., et al., 1998, ApJ 502,

L173

Walker C.K., Maloney PR., Serabyn E., 1994, ApJ 437, L127
Walsh A.J., Burton M.G., Hyland A.R., Robinson G., 1998, MNRAS
301, 640

Wilson T.L., Snyder L.E., Comoretto G., Jewell RR., Henkel C, 1996,
A&A 314, 909

Winnewisser G., Hocking W.H., Gerry M.C.L., 1976, J. Chem. Rhys.
Ref. Data 5, 79

Wright M., Sandell G., Wilner D.J., Plambeck R.L., 1992, 393, 225
Zinchenko 1., 1995, A&A 303, 554

Zinchenko 1., Mattila K., Toriseva M., 1995, A&AS 111, 95
Zinchenko I., Pirogov L., Toriseva M., 1998, A&AS 133, 337